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Abstract 

We use a modified SU(2) chiral sigma model to study nuclear matter at high density 
using mean field approach. We also study the phase transition of nuclear matter to 
quark matter in the interior of highly dense neutron stars. Stable solutions of Tolman- 
Oppenheimer-Volkoff equations representing hybrid stars are obtained with a maximum 
mass of 1.69M0, radii around 9.3 kms and a quark matter core constituting nearly 
55-85% of the star radii. 
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1 Introduction 



It has been argued that nuclear matter undergoes a phase transition to quark matter 
at high densities and/or high temperatures. It is expected that the high temperature 
limit has interesting consequences in heavy ion collision and/or in cosmology, but high 
baryon density behavior is important to the study of neutron stars . 

The quark structure of hadrons implies that at sufficiently large nuclear densities 
nuclear matter should convert into quark matter. The density at which transition occurs, 
is beheved to be a few times nuclear matter density . The lattice calculations indicate 
that for nonzero quark masses the phase transition may be weak first order or second 
order^. Most of the model calculations find it to be first order. Thus for large enough 
mass of neutron star, its core may consists of quark matter. In addition, if the phase 
transition is first order, a part of the core may consists of mixed phase of quark and 
nuclear matter. Kapusta and others^'^ have used a nonlinear Walecka model for the 
nuclear phase and a bag model for the quark phase . They have found that if the hadron- 
quark transition density ub > 4no(no is the nuclear matter density =0.153fm~^),then it 
is quite unlikely that stable stars with quark interior exist. Ellis et al.^ have also studied 
the possibility of a second order hadron-quark phase transition. In such cases, one 
needs to include an additional phenomenological parameter. A number of studies using 
different models have also been undertaken^"^ showing a first order phase transition. 

In all the above works, (baryon chemical potential)is the only conserved charge 
and pressure remains constant in the mixed phase. Glendenning et al.^ and Burgio 
et al.^ have considered both hb and he (electron chemical potential ) as conserved 
charges. This has the consequence that pressure varies continuously for all mixtures of 
the two components (hadron and quark) throughout the mixed phase. 

In the present work, we have used a modified SU(2) chiral sigma model(MCH)^° for 
hadronic matter since chiral model has been very successful, as such, in describing high 
density nuclear matter. The importance of chiral symmetry^^ in the study of nuclear 
matter was first emphasized by Lee and Wick^^ and has become over the years, one 
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of most useful tools to study high density nuclear matter at the microscopic level. The 
nonlinear terms in the chiral sigma model give rise to the three-body forces which become 
significant in the high density regime-*^^. Further, the energy per nucleon at saturation 
needed the introduction of isoscalar field^"^ in addition to the scalar field of pions^^. We 
also include the interaction due to isospin triplet p- vector meson to describe the neutron 
rich-matter^^. 

The modified SU(2) chiral sigma model^^ considered by us includes two extra higher 
order scalar field interaction terms which ensures an appropriate incompressibility of 
symmetric nuclear matter at saturation density. Further, the equation of state(EOS) 
derived from this model is compatible with that inferred from recent heavy-ion collision 
data^^. 

In our work, we consider the baryon chemical potential //b as the only conserved 
charge. Consequently pressure remains constant in the mixed phase region. A first 
order phase transition between beta stable nuclear matter and quark matter is indicated. 
Taking the existence of such a phase transition between nuclear matter and quark matter 
as a guide, we solve the Tolman-Oppenheimer-Volkoff(TOV) equations with appropriate 
nuclear matter^° and quark matter^^ equations of state and find the hybrid stars to 
consist of a quark-matter core with the nuclear matter forming the crust. 

This paper is organised as follows. In sec. 2, we present the equation of state for 
nuclear matter. The quark matter equation of state is discussed in sec. 3. In sec.4, we 
discuss the structure of hybrid star. We discuss and summarise our results in sec. 5. 

2 Nuclear matter Equation of State 

The modified SU(2) chiral sigma model considered by us is described by the Lagrangian 
density^°. 





\B 




6m^ 
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In the above Lagrangian, F^^ = d^uj^ - d^uj^, G^^ = d^p^ - d^p^ and x = (tt^ + (7^)^/^, -0 
is the nucleon isospin doublet, tt is the pseudoscalar-isovector pion field, a is the scalar 
field and o;^, is a dynamically generated isoscalar vector field, which couples to the 
conserved baryonic current = ip^jf^ip. pli is the isotriplet vector meson field with mass 
nip. B and C are constant coefficients associated with the higher order self- interactions 
of the scalar field . 

The masses of the nucleon, the scalar meson and the vector meson are respectively 
given by 

m = g^xo, rria = \/2Xxo, = g^xo (2) 

Here Xq is the vacuum expectation value of the a field , gi_^, gp and g^ are the coupling 
constants for the vector and scalar fields respectively and A = (rr4 — rn^)/{2f^), where 
m^r is the pion mass , /tt is the pion decay coupling constant . 

Mean-field approximation has been used extensively to obtain field theoretical equa- 
tion of state for high density matter^^. Using this approximation, the equation of motion 
for isoscalar vector field is 

^^0 = ^ (3 

and the equation of motion for the scalar field in terms of |/ = ^ is of the form^'^ 

where m* = ym is the effective mass of the nucleon and the coupling constants are 



mi mi 

The baryon number density ub — Up + Un — q;;^ , where /c/ is the Fermi momentum 
and 7 is the spin degeneracy factor which is equal to 4 and 2 for nuclear and neutron 
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matter respectively. The equation of motion for p , in the mean field approximation 
gives ■'^^ 

pI = {gp/2ml){np-nn) (6) 

At high densities the interior of neutron stars composed of asymmetric nuclear matter 
with an admixture of electrons. The concentrations of neutrons,protons and electrons 
can be determined using conditions of beta equilibrium and electrical charge neutrality^°. 

(here pi is the chemical potential of the particle species i) . In our analysis we find the 
ratio Up/un to he in the range 0.003-0.08 for ub taking value from 0.02 fm~^ to 0.25 
fm~^. We have included the interaction due to isospin triplet p- meson in Eqn-1 for 
describing neutron-rich matter. The symmetric energy coefficient that follows from the 
semi-empirical nuclear mass formula is^°'^^ 

where Cp = gl/ml . 

The diagonal components of the conserved total stress tensor corresponding to the 
Lagrangian(Eqn.l) together with the mean field equation of motion for the fermion field 
and a mean-field approximation for the meson fields is used to calculate the equation of 
state . The total energy density (e) and pressure(P), for the neutron rich nuclear matter 
in /3-equilibrium is given by^° 



^ mM-y^f ^ _ C _ 



n,p,e • 



■"^(^ - y'^' + - ^(1 - 



n,p,e- 

E 



The energy per nucleon is ^ = and the chemical potential is p — (P + e)/nB- 
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Figure 1: Pressure (P) vs. baryon number density(nB) for nuclear matter 

The values of five parameters C^, Cuj, Cp, B and C occurring in the above equa- 
tions are obtained by fitting with the saturation values of binding energy /nucleon (-16.3 
MeV), the saturation density (0.153 fm~^), the symmetric energy (32 MeV), the effec- 
tive(Landau) mass (0.85M)^^, and nuclear incompressibility (~300 MeV), in accordance 
with recent heavy-ion collision data^^ are — 1.999 fm^, = 6.8157 fm^, Cp — 4.661 
fm^, B = -99.985 and C = -132.2456. 

The pressure as a function of number density of nuclear matter is presented in Fig.[l]. 
The solid curve(MCH) corresponds to the model considered by us, whereas the dotted 
curve(CH) corresponding to the original chiral sigma model^^ is presented for compar- 
ison. The equation of state for the present model is found to be softer with respect to 
the original one. 
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3 Quark Matter Equation of State and Phase Tran- 
sition. 



Several authors^^"^^ have studied the possible existence of quark matter in the core of 
neutron stars/pulsars. Densities of these stars are expected to be high enough to force 
the hadron constituents or nucleons to strongly overlap thereby yielding quark matter. 
Since the distance involved is small, perturbative Quantum Chromodynamics (QCD) is 
used to derive quark matter equation of state . We consider here the quark matter EOS 
which includes u, d and s quark degrees of freedom^^'^^ in addition to electrons. We 
have taken the electron, up and down quark masses to be zero^^ and the strange quark 
mass is taken to be 180 MeV. In chemical equilibrium, one has iJ,d= l^s = l^u + A*e which 
can be written in terms of baryon and electric charge chemical potentials as^'^ 



1 1 



He = -He (8) 



„2 



The pressure contributed by the quarks is computed to order a — where g is the 
QCD coupling constant. The electron pressure is^^ 



(9) 



and the pressure for quark flavor /, with / = u, d or s is 2.18.25 



(10) 

Where kf — (//^ — m^)^/^ is the Fermi momentum. The total pressure, including the 
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bag constant B is given by 



.18 



(11) 



In the above equations ij,b and ij,e are only two independent chemical potentials, /xg 
is adjusted so that the matter is electrically neutral i.e dP/djiB — . The baryon 
number density {ub) and the energy density (e) for quark matter can be derived using 
the thermodynamic relations^^ 



Now we shall study the possible scenario of phase transition from nuclear matter 
to quark matter . Gibb's criteria is used to determine the phase boundary of the co- 
existence region between the nuclear and quark phase. The critical pressure and the 
critical chemical potentials are determined by the condition 



We have taken typical values ct^ = 0.5, 0.6 and the bag constant B — (155 MeV)^ and 
(150 MeV)^, which are reasonable values to calculate pressure in the quark sector^. 

We have plotted pressure versus chemical potential for beta stable nuclear matter 
and quark matter in Fig. [2]. The solid line is shown for the beta-stable nuclear matter 
with MCH parameter set. The dotted and dashed lines correspond to quark matter with 
as — 0.5 and ag — 0.6 with bag pressure B — (155 MeV)^ respectively . We find that 
there exist phase transition points for nuclear matter at different pressures and chemical 
potentials. The transition point {Pcrit., I^cru.) for ag — 0.6 is (101 MeV/fm^, 1237 MeV) 
and that for — 0.5 is (43 MeV/fm^, 1111 MeV). Fig. [3] shows the phase diagram with 
a different bag pressure B = (150Mey)^. Here the transition point {Pent., fJ-crit.) for 
= 0.6 is (27 MeV/fm3, 1064 MeV) and that for = 0.5 is (7 MeV/fm^, 993 MeV) . 
Thus, we find though the transition point increases with the increase of bag pressure 
and as, with one of them remaining constant, the dependence on ag is more sensitive. 

Now considering a typical transition point {Pent., fJ'crit.) for B=(150MeV)'^, as= 0.5, at 
the critical pressure, the energy densities for the quark matter (e^J^t.) and nuclear matter 



Ub — dP/d/iB and e 



Pnm{fJ'B) — PqmifJ'B) 



(12) 
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Figure 2: Pressure (P) vs. chemical potential for nuclear matter and quark matter with 
various ag at constant bag pressure (B) 

(CrTt.) sectors are found to be 387 MeV/fm^ and 252 MeV/fm^ respectively. The baryon 
number densities corresponding to the critical ixb in quark matter is = 0.37/m~^ 
and that in nuclear matter is = 0.27 fm~^ which is reasonably an order of magnitude 
(about 1.8 to 2.5 times) higher than the nuclear matter density. The discontinuity in the 
number density as well as energy density indicates a first order phase transition. This 
phase transition from nuclear matter to quark matter obviously implies that the interior 
of neutron star consists of quark matter. We investigate this possibility further in the 
next section. 
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4 Hybrid Stars. 

Having established the existence of a phase transition, we now proceed to study the 
structure of a hybrid star. For the description of neutron star which can generate 
curvature in the space time geometry due to high concentration of matter, one has to 
apply Einstein's general theory of relativity. The space-time geometry generated by a 
spherical neutron star, described by the Schwarzschild metric can be represented in the 
form^^ 

^ _e-«(it2 ^ [1 _ 2M{r)/r]-^dr^ + r^[dQ^ + sin^Qd(t)\ (13) 

The Tolman-Oppenheimer-Volkoff (TOV) equations which determine the star structure 
and the geometry, in dimensionless forms^^'^^, are given by 

dPjrro) ^ ^ [e(fro) + P{rro)][M{fro) + A7raf^P{fro)] 
df ~ f2[l - 2GM(rro)/f] 

M(fro) = 47ra / df'r' e(fVo) , (15) 
Jo 

and the metric function, z/(r), relating the element of time at r = oo is given by^^ 

du{frQ) _ ^^ [M(rro) + 4:7iaf^P{fro)] 
df ~ f2[l - 2G'M(fro)/f] 
The following substitutions have been made in above Eqns.(14-16). 

e = e/ec, P = P/ec, f = r/ro, M = M/Mq (17) 

Here , 

a^e.roVMe, (18) 
with /i = 197.327 MeV fm and Tq = 3x 10^^ fm. 

The quantities with hats are dimensionless in above equations. The gravitational con- 
stant G = 6.7079 X 10-^^ MeV-2. 

For complete calculation of a stellar model, one has to integrate Eqns.(14-16) from 
the star's center at r = with a given central density Ec as input until the pressure P(r) 
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at the surface vanishes. With any reasonable central energy density, we expect that at 
the center we shall have only quark matter. Hence we shall be using here the equation 
of state for quark matter through Eqn.(lO) with P(0) = P(ec)- We then integrate the 
TOV equations until the pressure and density decrease to their critical values at r = Tc- 
For r > Tc, we shall have equation of state for /^-stable nuclear matter where pressure 
will change continuously but the energy density will have a discontinuity at r = Tc . The 
TOV equations with nuclear matter equation of state are continued until the pressure 
vanishes which defines the surface of the star. This completes the calculations for stellar 
model for a hybrid neutron star, whose mass and radius can be calculated for different 
central densities. 

Fig. [4] shows behavior of pressure versus number density in the vicinity of a first order 
phase transition in a system having one chemical potential corresponding to conserved 
baryon number^. Points labeled H and Q mark the end of the hadronic and beginning 
of the quark phase, the intervening region representing the mixed phase. The two equal 
pressure points at the opposite ends of the mixed phase are mapped onto the same radial 
point in the star. These aspects are also illustrated in Fig. [5]. 

Through out the calculation of hybrid star we have used ag — 0.5 and B = (ISOMeV)^ 
for quark sector. The energy density profile obtained from Eqns. (14-16) are plotted in 
Fig. [5] for central energy density ec— 450 MeV/fm^ and ec = 1700 MeV/fm^. For core 
energy densities greater than the critical energy density (~ 387 MeV/fm^) the core 
consists of quark matter. As we go away from the core towards the surface through 
TOV equations, when the critical pressure is reached, the density drops discontinuously 
indicating a first order phase transition. Thus for central density of 450 MeV/fm^ such 
a star has a quark matter core of radius 5.5 kms with nuclear matter crust of about 2.2 
kms , whereas for ec = 1700 MeV/fm^, the quark matter core radius is 8.2 kms with 
nuclear matter crust of about 0.6 kms. Hence it is clear that if we take smaller central 
energy densities then nuclear matter is expected to be more abundant in a hybrid star . 

We have plotted in Fig. [6] the mass of hybrid star as a function of central energy 
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density to examine the stability of such a star. Hence we have two branches of solutions. 
Pure neutron star at lower densities Ec < and hybrid stars at central densities 
Ec > CgJ^. Taking into account the stability of such stars under density fluctuations 
require dM/dCc > 0^^. As may be seen from the figure, dM/de^ becomes negative 
around 1700 MeV/fm^ after which it may collapse into black holes^^'^^. This yields the 
maximum mass of hybrid star as M ~ 1.69 Mq. Fig. [7] shows the mass as a function 
of radius obtained for different central densities varying in the range of 450 MeV/fm^ 
to 1700 MeV/fm^ for such a star which indicates the maximum radius to be around 9.3 
kms. To check the stability of our result against the errors in values of symmetric energy 
(32±6 MeV) and the nuclear incompressibility (300±50 MeV), we have also computed 
the coupling constants (Cp, C^^, Cu,B and C) for four cases of symmetric energy value 
32 MeV with incompressibihty 250 and 350 MeV and the incompressibility of 300 MeV 
with symmetric values 26 and 38 MeV. The maximum values for M and R are found 
to be M = 1.69to;ooiMQ and R = 9.31q;o5 kms. Such small effect of symmetric energy 
and incompressibility on maximum values of M and R can be attributed to the fact that 
only about 15% of the radius of the hybrid star contains nuclear matter. 

We also calculate the surface gravitational redshift Zg of photons which is given by^^ 

1 (19) 



'(1 - 2GM/R) 

In Fig. [8] we have plotted Zg as a function of M/Mq . In this context it may be 
mentioned here that our result for the surface redshifts lie in the range of 0.2 to 0.5 as 
determined from gamma ray bursters^°. 

We then compute the relativistic Keplerian angular velocity fl^ given by^^ 



= 0.72, 



M/Mq 



lO^sec-i ^"^\\{R/10kmy 
as for neutron stars. Fig. [9] shows our result for variation of relativistic Keplerian 

angular velocity as a function of M/Mq for such a star. We find that the flk has an 

inverse relationship after an initial increase with the mass of the star beyond 1.69M0^. 

11 



This indicates that we can not have mass of hybrid star more than about 1. 69 Mq. We 
also observe that the maximum value of Jlj. being near about lO^sec"^ implies the time 
period to be less than 0.5 millisecond. 

5 Conclusions 

Within the formalism of a modified SU(2) chiral sigma model (Eqn.l) we found that a 
first order phase transition exists between the nuclear phase and quark phase at density 
of about two to three times the nuclear matter density. Quark matter equation of state 
has the parameters: m^, ag and B. Our results show that the critical parameters while 
increasing with increasing value of cts and B, exhibit greater sensitivity with respect to 

The phase transition from nuclear matter to quark matter indicates that the core of 
a neutron star consists of quark matter. To obtain a stable hybrid star solution, we have 
solved TOV equations using appropriate equations of state with a given central energy 
density Cc and have taken ag — 0.5 and B = (150Mel^)^ for quark matter EOS. It is 
observed that a stable hybrid star with a quark core and a nuclear matter crust exist 
upto ec ~ 1700 MeV/fm^ beyond which instabihty is indicated. For ec varying from 450 
MeV/fm^ to 1700 MeV/fm^ the mass of hybrid star varies from 0.5 to 1.69 Mq and the 
radius from 7.7 to 8.8 kms with a quark core of about 5 to 8 kms respectively. Our 
results, thus, indicate that the bulk of hybrid star is composed of quark matter with a 
crust of nuclear matter. We find that the maximum mass and radius of hybrid stars to 
be about 1.69 Mq and 9.3 kms respectively depending on the values of the parameters 
used in the model. If we compare these mass and radius values of the hybrid stars with 
those obtained using this modeP° for pure neutron star (such as M = 2.1Mq and R — 
12.1 kms), we find that the hybrid stars are more compact than a normal neutron star. 
This result is in exact agreement with the general expectation since the EOS of quark 
matter forming the core of hybrid star is supposed to be softer than that of neutron 
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matter because of QCD asymptotic freedom. The greater compactness of the star also 
leads to smaller time period lying in the submillisecond range as obtained by us. It is 
also observed that the surface garvitational redshift and relativistic Keplerian angular 
velocity of the hybrid stars can not increase beyond M/Mq = 1.69 ; showing a decrease 
with increase in M/Mq beyond this value. 
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Figure 3: Pressure (P) vs. baryon chemical potential for nuclear matter and quark 
matter with various with different bag pressure (B) 
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Figure 4: Pressure (P) vs. number density(nB) of hybrid star. 



17 




R (km) 

Figure 5: Energy density profile for central densities €c — 1700 MeV/fm^ (solid line) 
andec = 450 MeV/fm^ (dashed line). 
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Figure 6: The mass(M/MQ) of the hybrid star as a function of central energy density(ec). 
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Figure 7: Mass as a function of radius for the hybrid star. 
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Figure 8: The surface gravitational redshift(Zs) as a function of star mass(M/MQ) 
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Figure 9: The Keplerian angular velocity(Qfe) as a function of star mass(M/MQ) 
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